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Macroscopic entanglement witnesses have been put forward recently to reveal nonlocal quantum 
correlations between individual constituents of the solid at nonzero temperatures. Here we apply a 
recently proposed universal entanglement witness, the magnetic susceptibility [New J. Phys. 7, 258 
(2005)] for the estimation of the critical temperature T c in the nanotubular system Na2Vs07 below 
which thermal entanglement is present. As a result of an analysis based on the experimental data 
for dc-magnetic susceptibility, we show that T c m 365 K, which is approximately three times higher 
than the critical temperature corresponding to the bipartite entanglement. 
PACS numbers:75.30.Et, 75.10.Pq, 03.67.Mn 
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I. Introduction 

Entanglement lies in the heart of quantum mechanics, 
revealing the existence of nonlocal correlations that can- 
not be described by classical physics alone [lj. Because 
of such strong correlations over and above the classical 
one, entangled states have been recognized as a crucial re- 
source for quantum information processing 2j. Thus for, 
experimental demonstrations of entanglement have been 
confined yet mostly to atomic scales. This is due to the 
subtle nature of quantum entanglement, since extending 
the size of the physical object increases the number of 
constituents, and the number of degrees of freedom that 
can interact with the environment, which is responsible 
for reducing and eventually destroying entanglement in 
the system. 

However, it was demonstrated recently that entan- 
glement can affect the macroscopic properties (such as 
heat capacity or magnetic susceptibility) of an insulat- 
ing magnetic compound LiHo x Yi_ x F4 jjj, but at a very 
low (below 1 K) temperature. The basic reason for it is 
that magnetic susceptibility is proportional to the two- 
site correlations between the spins in the magnetic solid, 
which can be higher than the one allowed by classical 
physics alone, when entanglement is present in the ma- 
terial as well. In the light of the result of Ghosh et al 
that entanglement can have a significant effect in the 
macroscopic world, the question arises whether we can 
experience entanglement at higher temperatures as well. 

As a next step, Brukner et al showed that for un- 
coupled dimers at zero magnetic field the magnetic sus- 
ceptibility serves as an entanglement witness (EW), i.e., 
an observable which can distinguish between entangled 
and separable states in a quantum system |5j . It has been 
demonstrated |j] that the experimental data on magnetic 
susceptibility in the spin-1/2 alternating bond antiferro- 
magnet CN (Cu(N03) 2 2.5D20) Q implies the presence 
of entanglement below the temperature 5 Kelvin. 

On the other hand, Bose and Tribedi showed, based 
on energy considerations 8], that the magnetic suscep- 
tibility as an EW can be applied to compounds which 
can be modeled as a collection of independent clusters 



of N spins, when the Hamiltonian is isotropic in spin 
space. As a special case, Bose and Tribedi considered the 
5 = 1/2 polyoxovanadate AFM compound V12, which 
material can be treated as a collection of independent 
5 = 1/2 tetramers and estimated that the critical tem- 
perature below which entanglement exists in the system 
is T c w 25.4 K. 

Other thermodynamical observables, such as internal 
energy, magnetization, or heat capacity |sl. 1^. ITol ITU IT2I 
Il3l ] , have also been proposed for EW but either they rely 
on a model Hamiltonian of the system (and usually not a 
directly measurable quantity) or are applicable only for 
relatively low temperatures |13| . 

However, recently Wiesniak et al 01 constructed an 
EW based on the sum of magnetic susceptibilities mea- 
sured along the three orthogonal spatial directions. This 
entanglement witness is applicable to systems without 
full knowledge of the specific model of the solid. The 
separability criterion of Wiesniak et al 0] for N spin-l 
particles reads as follows: 



X = Xx 



Xy 



Xz>7f l i 



(1) 



with C = (gfJ-B) N/k, where g is the g-factor, [is is 
the Bohr magneton, k is the Boltzmann constant, and 
N denotes the number of particles in the volume. Fur- 
ther in the right-hand side of JIJ, x denotes the sum of 
magnetic susceptibilities per volume measured along the 
orthogonal x, y, and z directions. In the isotropic case 
for spin-1/2 particles inequality simplifies to the in- 
equality 



Xz > 



CI 
T6 



(2) 



Therefore, if Xz < ^\ the isotropic solid state system 
contains entanglement between the individual spin-1/2 
systems. Note, that Xz entanglement witness in (0 is the 
same as the one in Ref. Q, although that has been in- 
troduced only for uncoupled dimers using a simpler argu- 
mentation. Since the entanglement witness x of Wiesniak 
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et al is completely general, the entanglement crite- 
rion (J2J can be applied for a large class of quantum spin 
systems, as it neither rely on the existence of small spin 
clusters or on the full knowledge of the system's Hamil- 
tonian. 

The low-dimensional quantum magnets have recently 
been in the focus of both experimental and theoreti- 
cal research. A large class of these materials are the 
spin-1/2 quantum magnets with chain, ladder, or planar 
geometries. In this respect new possibilities emerged, 
when Millet et al 0] recognized that the variation of 
some two-dimensional compounds results in a nanotubu- 
lar structure, and as a result the first transition-metal- 
oxide-based nanotubular system Na2Vs07 was synthe- 
sized. In this system the VO5 pyramids share edges and 
corners to form a nanotubular structure with Na atoms 
located inside each individual tube. The ab initio analy- 
sis in Ref. 16] showed that Na2Vs07 can be described as 
formed by nanotubes consisting of weakly coupled nine- 
site rings. 

The main purpose of the present paper is to demon- 
strate that the entanglement witness \z m criterion (0 
implies entangled states in the system Na2Va07 at the 
room temperature range, and we estimate the critical 
temperature T^ xp m 365 K. The validity of this result 
is also supported by additional calculations. In par- 
ticular we find that neglecting the inter-ring couplings 
as a first approximation, the inequality |(SjSj)| < 1/4, 
where i,j denotes spin sites within the ring, may indi- 
cate pairwise entanglement. Based on the parameters of 
the model Hamiltonian of the system Na2Vs07 which 
are derived from ab initio molecular structure calcula- 
tions fl6| . we calculate the critical temperature corre- 
sponding to |(SiSj)| = 1/4, and obtain Tf alT « 123 K, 
below which pairwise entanglement is present in the sys- 
tem. The three times smaller temperature bound for 
the bipartite entanglement indicates that a large portion 
of entanglement is stored in correlations between more 
than two spins within the nine-site ring. Therefore we 
conclude that the entanglement above the temperature 
rppair j g multipartite in nature. To account for this effect 
we analyze the symmetric trimer (described by S = 1/2 
Heisenberg antiferromagnetic (AFM) isotropic Hamilto- 
nian), since by small temperatures the nine-spin cluster 
model of the Na2Va07 system corresponds to a three- 
site ring. Although the symmetric trimer is known not 
to have pairwise entanglement both at T = and at 
finite temperatures |l7j . we find multipartite entangle- 
ment witnessed by criterion J5J in the system. We ar- 
gue that the "missing entanglement" in the case of the 
symmetric trimer, not explored by the |(S;Sj}| < 1/4 
condition, appears in the nanotubular system Na2Va07 
as well, and the existence of this multipartite entangle- 
ment causes the large difference in the values of critical 
temperatures T^ xp and Tf alT ' . Further, we claim that 
due to the weakly coupled nine-site rings with respect 
to the intraring couplings, the entangled state itself in 
the nanotubular system is shared with a good approxi- 



mation among the nine spins within a ring. Therefore 
the entangled state is confined to the rings of the nan- 
otube, and does not extend to macroscopic sizes. This 
fact may partially explain why the entanglement in this 
system persists up to such high temperatures. 

II. Entanglement witnesses based on two-site spin 
correlation functions 

When a system is in its thermal equilibrium under a cer- 
tain temperature T, it can be described by the density 
operator p (T) = Z^ 1 cxp (—f3H), where H is the Hamil- 
tonian, Z is the partition function, and /3 = 1/kT. A 
thermal state remains entangled up to a critical temper- 
ature T c above which the state becomes separable, i.e., 
the amount of entanglement vanishes. 

In order to measure the thermal entanglement based 
on the thermal density matrix we need to know 
both the eigenvalues and eigenvectors of the Hamilto- 
nian. On the other hand, bulk properties of solids 

such as internal energy (u = — Z^ 1 -^-^, magnetiza- 
tion ((M z ) = — {Z(3)~ l -7^^j , or magnetic susceptibility 

^Xz = 9 qb* ' J 1 can be derived from the partition func- 
tion alone; thus these quantities are defined merely by 
the eigenvalues of the system H . 

The latter relationship between magnetization and 
magnetic susceptibility has given rise to the notion of 
thermodynamical complementarity [Tsj . On the other 
hand, if [M, H] = the equality 

Xz = = Jl {(Ml) - (M z f) (3) 

holds linking the fluctuation of magnetization to the 
macroscopic observable \z- Here A 2 M Z denotes the vari- 
ance of the magnetization, while the constant C has the 
same meaning as in inequality Q and N is the number 
of particles in the volume. On the other hand, noticing 

N 

that M z — Sf , formula J2Jl also connects the mag- 

i=l 

netic susceptibility to its microscopic origin in the form 
of spin-spin correlation functions. 

Now we represent two entanglement witnesses from 
the literature, based on spin-spin correlation functions, 
which are crucial to construct macroscopic entanglement 
witness from the magnetic susceptibility through equa- 
tion This promise to reveal thermal entanglement in 
macroscopic systems even when the eigenvectors of the 
system's Hamiltonian are not known. For any separable 
states of N spins of length I, that is, for any convex sum 
of N product states, p = Yl ^kp\®p\ ® • • • <8> , one has 

, . k 

El 

N 

El(SA + i)|<iVi 2 , (4) 

i=l 

which especially for two qubits (at sites i and j) yields 
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the inequality 



J(SiS 2 + S 2 S 3 + S3S1), with J > thus criterion 



(5) 



The other inequality for separable states of N spins of 
length I is based on entan glem ent detection applying the 
sum uncertainty relations [TiL , 



A 2 M X + A 2 My + A 2 M Z > Nl . 



(6) 



Now we restrict ourselves to Heisenberg spin-1/2 lattices 
at zero field with isotropic, but in general inhomogeneu- 
ous coupling constants Jy , 



H 



(7) 



Due to the isotropy of the Hamiltonian in spin space, the 
magnetization vanishes in all three orthogonal directions, 
and thus A 2 M a = (M 2 ) for a — x,y, z. Using the rela- 

N 

tion M a — Sf and assuming spin-1/2 lattices Eq. © 



simplifies to 



JV 



> 



N 
2~ ' 



(8) 



concur- 



From now on for brevity we denote the entanglement 
witness corresponding to criterion (JSj by EW1, and by 
EW2 of the criterion © . Note that combining inequal- 
ity © with © implies the criterion (see Wiesniak et 
al 14] ). On the other hand for the isotropic Heisenberg 
Hamiltonian 0) magnetization M commutes with the 
Hamiltonian H , [M, H] = 0, thus the witness Xz in crite- 
rion J2J is valid. As a consequence the power of EW2 in 
(|5|) to signal entanglement is equivalent with the witness 
Xz in criterion |[2J. 

Wang and Zanardi |2(| showed that for the Heisen- 
berg Hamiltonian Q the concurrence |2l|, which 
gives a measure of pairwise entanglement between 
qubits at site i and site j, becomes C, 
1/2 max [0,2|<S?S?)| - (SfSj) - 1/4]. The 
rence CV,- = corresponds to an unentangled state. Re- 
garding to the isotropy of the model the a — x, y and z 
components of the spin-spin correlation functions are all 
equal, and assuming that (SjSj) is negative, we arrive 
at dj = 1/2 max [0, - (SiSj) - 1/4]. Thus, whenever 
(SjSj) < holds, pairwise entanglement between qubits 
i and j vanishes for |(SjSj)| < 1/4. Since the condition 
above is just the separability criterion (J5J, this implies 
that EW1 detects pairwise entanglement between any 
two qubits of a system with Hamiltonian when the 
corresponding spin-spin correlation has negative value. 

Now we compare the power of EW1 and EW2 de- 
tecting entanglement in a thrcc-qubit ring interacting 
via the AFM isotropic Heisenberg Hamiltonian and 
show that EW2 can be more effective. For the AFM 
symmetric trimer the Hamiltonian is given by H = 



gives £ (SiSj) = 3 ((Sf > + 2 (S,S i+1 )) > § and conse- 

quently (SjSj+i) > —1/8, for which the AFM interaction 
implies the relation |(SjSi+i)| < 1/8 for separable states. 
Thus in this case the condition |(SiS i+ i}| < 1/8 of EW2 
obviously implies a lower bound on the separability than 
the condition |(SjS i+ i)| < 1/4 of EW1. Since direct cal- 
culation shows that for the ground state of the symmet- 
ric trimer (SjSj+i) = —1/4, EW2 indicates entangled 
ground state, whereas EW1 could not reveal entangle- 
ment. The result that there exists entanglement in the 
symmetric trimer over the two-site entanglement is sup- 
ported by the definition of 1-tangle t\ as well 22] , which 
quantifies the entanglement of one site with the rest of 
the chain at T — 0, and is given as t\ — 4det (pi), where 

Pl = (i + £ {Si) f 7 "^ /2 is the one-site reduced density 

matrix, and a a are the Pauli matrices with a = x, y, z. 
In terms of spin expectation values {Sf), t% takes the 
form T\ = 1 — 4£ (Sf) 2 . The vanishing of n implies 

a 

that there is no entanglement in the ground state, i.e. 
the state is factorized. The isotropy of the symmetric 
trimer (Sf) = {Sf) = {Sf) = entails that n = 1, 
which means that indeed the ground state of the sym- 
metric trimer is entangled. Although in the above anal- 
ysis we treated merely the three qubit isotropic Heisen- 
berg ring, we conjecture that the results obtained above 
can be generalized to systems with general Heisenberg 
Hamiltonian Q for cluster size TV > 3 and for different 
but positive Jij values as well. We intend to support this 
conjecture in the next sections with the analysis of the 
quantum magnet Na2Vs07, which can be considered as 
a nine-leg nanotubular system, and is well described by 
the Hamiltonian (JJJ with J,j > exchange constants. 

III. Magnetic susceptibility of the quantum mag- 
net Na 2 V 3 7 

Na2V3C>7 is a material whose structure may be consid- 
ered as composed by nanotubes oriented along the c 
axis of the crystal lattice. A detailed description of the 
structure of this compound can be found in Ref. |15| . 
Following the discussion and ab initio results given in 
Ref. [16| we construct the model Hamiltonian of the sys- 
tem Na2V3C>7, which is given by 



( 9 ^ 



(9) 



where periodic boundary conditions were imposed 
(Si + g = Si), and the weak inter-ring couplings as a first 
approximation were neglected in the model. Thus in this 
model we treat the system as the collection of indepen- 
dent nine-site rings. The optimal parameters, neglecting 
also the differences in between the nearest neighbor (NN) 
and in between the next nearest neighbor (NNN) intrar- 
ing interactions, are given by J\ = 200 K and Jo = 140 K 
[la ]. Due to the structure of the material the ring 
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FIG. 1: Illustration of exchange couplings Ji, J2 in the spin- 
1/2 nine-leg ring model of the nanotube Na2Vs07. 



is partially frustrated, where J\ = for i — 1,4,7 and 
J% = J2 for the indices i = 2, 3, 5, 6, 8, 9. The two kinds of 
couplings in the frustrated model are depicted in Fig. 1. 

The ground state of the model is fourfold degenerate, 
composed of two doublets, corresponding to the right and 
left chirality |23j . Since this model is invariant under 
translations i — > i + 3 it corresponds by low energies to 
a three-site ring. In Fig. 2 we display the experimentally 
observed inverse magnetic susceptibility of the Na2Va07 
at temperatures between 1.9 K and 315 K (data is taken 
from Ref . [24| ) . As a comparison we also reproduced and 
plotted in Fig. 2 the calculated inverse magnetic suscep- 
tibility for the partially frustrated model from Ref. [l(| 
with couplings J\ = 200 K and J2 = 140 K in the temper- 
ature range T = — 600 K. The temperature dependence 
of the magnetic susceptibility was obtained in Ref. 0] 
from the thermal state p(T) by exact diagonalization of 
the Hamitonian @ • The nice agreement between the ex- 
perimental and calculated values supports the validity of 
the ab initio derived parameters J\ and J2 in the model. 

IV. Calculation of thermal entanglement in the 
system Na2Va07 

In this section we give a lower bound on the critical tem- 
perature below which the thermal state of the system 
Na2V3C>7 is entangled. The two estimations for this tem- 
perature limit will be calculated based on entanglement 
witnesses EW1 and EW2 (in criterions (JSJ and JSJl, re- 
spectively). EW1 is able to reveal pairwise entanglement 
between individual spins, while EW2 may reveal resid- 
ual entanglement which is not detected by EW1, as it 
was demonstrated in the example of symmetric trimer. 
The critical temperature related to the pairwise entangle- 
ment and detected by EW1 will be called T pa " ', while the 
possibly higher critical temperature indicated by EW2 is 
denoted by T^ xp . The relation between magnetic sus- 
ceptibility and two-site spin correlation functions links 
EW2 to the entanglement witness \z expressed by con- 
dition J2J). Therefore the critical temperature revealed by 
the witness \ z also applies to EW2, and in the following 
it is jointly denoted by T^ xp . 

Since the magnetic susceptibility of the system 
Na2V3C>7 is available from Sec. Ill, we next calculate 
the critical temperature T£ xp signaled by % z . Applica- 



o 



600 
500 
400 
300 
200 
100 




• experiment 




model 








f Entanglement 





100 



200 



300 l ° 

T(K) 



-365 K 



500 



600 



FIG. 2: Temperature dependence of the inverse magnetic sus- 
ceptibility in units of mol/emu of the system Na2Vs07. The 
filled dots correspond to the experimental data from Ref. [2^ |. 
The solid curve is the calculated curve by reproducing the 
results of Ref. [l6| , and the dotted line represents the macro- 
scopic entanglement witness Xz • The vertical solid line shows 
the critical temperature T^ xp ~ 365 K, i.e., the intersection 
point of the dotted and solid curves, below which entangle- 
ment exists in the nanotubular system Na2Vs07 . 



tion of criterion J2J) implies that entanglement is present 
in the system, when the inverse magnetic susceptibility 
per volume satisfies the inequality 



X: 



1 >6 



T 

C ' 



(10) 



where the concrete value of C = (g/is) N/k for the ma- 
terial Na2V3C>7 is given by C as 5.58 emu K/mol. The 
equality in (|10f> is represented by the dotted line in Fig. 2. 
The intersection point of the solid curve (calculated in- 
verse susceptibility versus T) and the dotted line provides 
the estimate of the critical temperature T c cxp « 365 K. 
One may note that experimental data was measured up 
to 315 K and therefore the estimation was obtained by 
the calculated curve, which is not so reliable. However, 
data of inverse susceptibility is larger than the right-hand 
side of O up to T = 315 K, thus T^ xp is definitely 
higher than 315 K, and therefore the experimental data 
in this temperature range cannot be described fully with- 
out taking into account entanglement. 

Now relying on the assumption that the partially frus- 
trated model 10 describes well the real physical system, 
we calculate a lower bound on the presence of entan- 
glement, witnessed by EW1 of criterion (J3J. If for any 
two sites of the ring of nine qubits the inequality 

|(SjSj)| > 1/4 holds, then the state of the pair (i,j), 
described by the reduced density matrix of those two, is 
entangled. Furthermore, regarding the result of Ref. |20|. 



5 



v 

0.125 
0.0 





Entanglement 

/ 




\j=2 


sS — — ■ 




\ 

\ 

/ 'v. 
- / x - 
/■ 









TP c' m ~i23K 200 300 4 oo 500 600 

T(K) 



FIG. 3: Temperature dependence of the absolute value of 
the two-site spin correlation function |(SiSj)|, for neighbor- 
ing (j = 2, 9) and second neighboring (j — 8, 3) sites. The 
solid curves correspond to nearest neighboring sites, while the 
chain curves represent spin-spin correlations between second 
neighbors. The horizontal dotted line represents the value of 
1/4 for the EW1. The intersection point of this line with the 
solid curve corresponding to j — 2 gives the critical tempera- 
ture T^ a%r » 123 K. The critical temperature is designated by 
a vertical solid line which defines the temperature range (left 
from the solid line) with pairwise entanglement in the system 
N a2 V 3 7 . 



if (SiSj) is negative we can state that the quantum cor- 
relation between the zth and jth sites is due to pair- 
wise entanglement. Fig. 3 shows the |(SiSj)| correlation 
functions versus T calculated in the temperature range 
T = - 600 K where the labels j = 2,9 and j = 8,3 
correspond to the first and second neighbors of site 1, 
respectively (see Fig. 1). We do not show spin-spin cor- 
relation functions for the third and fourth neighboring 
sites because of their small values. By the calculation 
of the correlation functions we used exactly the same 
model parameters as by the calculation of the inverse 
magnetic susceptibility in Sec. Ill so that the results in 
Fig. 2 and Fig. 3 would be consistent. Due to the sym- 
metry of the model system, there are two inequivalent 
neighboring (j = 2,9) and second neighboring (j = 8,3) 
sites. Numerically we obtained that apart from the case 
(S1S3) (i.e. j — 3), all the other correlation functions 
(i.e. j = 2, 9, 8) are negative in the whole range of 
T = — 600 K. On the other hand, it can be noticed 
in Fig. 3 that the only curve which intersects the value 
of 1/4, represented by the horizontal dotted line, is the 
function {S1S2) | - Since (S1S2) in the whole temperature 
range is negative, it implies that the estimated value of 
the critical temperature T£ aiT k 123 K refers to bipar- 
tite entanglement. According to the above analysis we 
conclude that EW1 is able to detect pairwise entangle- 
ment between six first neighboring sites, namely between 
(i,i) = (l,2), (2,3), (4,5), (5,6), (7,8), (8,9). 



We have also calculated the spin-spin correlation for 
the T = ensemble, where the ground state density ma- 
trix was assumed to be an equal mixture of all eigen- 
states, corresponding to the fourfold degenerate ground 
state. For the next nearest neighboring (NNN) spin pairs 
the values (SiS 8 ) = -0.011 and (S1.S3) = +0.1807 were 
obtained. As it can be observed in Fig. 3, these are just 
the limiting values of these NNN spin-spin correlation 
functions by T = 0, which indicates smooth spin-spin cor- 
relation functions near T = 0. The large difference in the 
values of (SiSg) and (Si S3) implies strong partial frus- 
tration in the system due to the large coupling J2 , which 
diminishes the spin-spin correlation between sites 1 and 
8 at T = 0. Although this effect may reduce the criti- 
cal temperature pertaining to the pairwise entanglement, 
qualitatively the result obtained in this section for the 
Na2V3C>7 is consistent with the result for the symmetric 
trimer we arrived at Sec. II. In the latter case the bound 
revealed by EW2 (|(SiS 2 )| = 1/8) was half of the bound 
corresponding to EW1 (|(SiS 2 )| = 1/4). On the other- 
hand, we can calculate the power of EW2 with respect to 
EW1 in the system Na2Vs07 as well. By eye inspection 
of the spin-spin correlation functions in Fig. 3, the value 
I {S1S2) I by the critical temperature T^ xp — 365 K ap- 
proximately yields | (S1S2) | = 1/10. This demonstrates 
the fact that despite of the frustration, there exists en- 
tanglement in the system Na2Va07 beyond the pairwise 
one. The residual thermal entanglement, stored in three- 
spin, four-spin, and higher-spin correlations, could not be 
detected by the |(SiS 2 )| = 1/4 bound of EW1, but was 
revealed by the macroscopic witness Xz (based on EW2). 

V. Conclusion 

In conclusion we have shown that the nanotubular system 
Na2Va07 is an ideal candidate to demonstrate thermal 
entanglement at the range of room temperature. The 
magnetic susceptibility as a macroscopic entanglement 
witness, proposed in Ref. 01 (and based on EW2 in 
criterion (JSJ) gives the T c exp 365 K lower bound on 
the critical temperature beyond which entanglement dis- 
appears, while the calculation of EW1 (of criterion (J5J) 
for the model Hamiltonian of the Na2Va07 system re- 
vealed pairwise entanglement up to the temperature 
repair ^ 223 K. This result indicates that pairwise en- 
tanglement is not the sole entanglement in the system, 
and the magnetic susceptibility \ z is able to exhaust a 
probably large portion of residual part of the total en- 
tanglement. However we have to mention that the total 
entanglement in this particular material does not extend 
to the whole system, but rather is confined to the rings 
of the tube. Although macroscopic observables, such 
as magnetic susceptibility, can be explained only on the 
ground of entangled states, the entangled state itself is 
not macroscopic but shared among nine spins within a 
cluster. We emphasize that while the calculation of the 
temperature bound T pair for the pairwise entanglement 
depends on the validity of the model of the real physi- 
cal system, the critical temperature T^ xp revealed by the 
entanglement witness \z does not require any assump- 
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tion about the explicit values of the exchange-coupling 
parameters of the system. Moreover, it does not rely 
on the existence of independent or weakly coupled spin 
clusters in the material, thus the result for T£ xp is reli- 
able even if our model assumptions do not hold exactly. 
The existence of large exchange coupling parameters are 
of great importance in order to demonstrate quantum 
entanglement at high temperatures in the solid. Since 
this requirement is fulfilled for a large class of materials, 



such as some low-dimensional [25| or finite [2(j quantum 
spin systems, we believe that entanglement as a crucial 
resource for quantum information processing in various 
types of solids might be harnessed at room temperature 
range as well. 
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